Oscillatory and monotonic modes of longwave Marangoni convection in a thin film 
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We study longwave Marangoni convection in a layer heated from below. Using the scaling k = 
0(\/Bi), where k is the wavenumber and Bi is the Biot number, we derive a set of amplitude 
equations. Analysis of this set shows presence of monotonic and oscillatory modes of instability. 
Oscillatory mode has not been previously found for such direction of heating. Studies of weakly 
nonlinear dynamics demonstrate that stable steady and oscillatory patterns can be found near the 
stability threshold. 
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Introduction. Marangoni convection in a liquid layer 
with upper free boundary is a classical problem in the 
dynamics of thin films and in the pattern formation [l], Q ■ 
In the pioneer theoretical paper, Pearson |3[ analyzed the 
linear stability of the layer with a nondeformable free sur- 
face. He considered two cases of thermal boundary con- 
ditions at the substrate: the ideal and poor heat conduc- 
tivity, when either the temperature or the heat flux are 
specified. In the latter case he found a monotonic long- 
wave instability mode for heating from below and zero 
Biot number Bi. For Bi 1 the critical wavenumber k 
is proportional to Bi 1 / 4 [lj. Many authors extended the 
analysis in order to include the deformation of the free 
surface. Review of analytical and numerical works can 
be found in [lj. In particular, several oscillatory modes 
were revealed; these modes were reported only for heat- 
ing from above. 

In the case of heating from below, a nonlinear analysis 
for ideally conductive substrate was performed in Ref. [J] : 
it was shown that the subcritical bifurcation occurs and 
instability with necessity results in film rupture. The be- 
havior of perturbations near the stability threshold was 
studied in [5| for the case of a poorly conductive sub- 
strate. Under assumption of large gravity, and, hence, 
small surface deflection, the amplitude equation was de- 
rived and the subcritical bifurcation was found. 

In this paper we demonstrate the existence of a new os- 
cillatory mode of longwave instability for the film heated 
from below. Using the scaling k — O(yTli), which was 
first suggested in Ref. [6], we derive a set of amplitude 
equations. Linear stability analysis gives both the mono- 
tonic and the oscillatory modes. Pattern selection near 
the stability threshold clearly demonstrates that insta- 
bility does not necessarily lead to rupture and that both 
steady and oscillatory regimes can be found experimen- 
tally within certain domains of parameters. 

Problem formulation. We consider a three-dimensional 
thin liquid film of the unperturbed height Ho on a pla- 
nar horizontal substrate heated from below. The heat 
conductivity of the solid is assumed small in compari- 
son with the one of the liquid, thus the constant vertical 
temperature gradient —A is prescribed at the substrate. 



(The Cartesian reference frame is chosen such that the x 
and y axes are in the substrate plane and the z axis is 
normal to the substrate.) 

The dimensionless boundary- value problem governing 
the fluid dynamics reads: 



Fv-Vv) = -Vp + V 2 v-Ge 2 , (la) 

T t + v • VT = V 2 T, V ■ v = 0, (lb) 
v = 0, T z = -1 at z = 0, (2a) 
S • n = (p - CaX) n - MV T (T\ z=h ) , V„T = -Bi T, 
h t = iv — v • V/i at z = h(x, y, t). (2b) 

Here, v = (u, w) is the fluid velocity (where u is the ve- 
locity in the substrate plane and w is the z-component), 
T is the temperature, p is the pressure in the liquid, £ 
is the viscous stress tensor, h is the dimensionless height 
of the film, e z is the unit vector directed along the z 
axis, n and r are the normal and tangent unit vectors to 
the free surface, respectively, K is the mean curvature of 
the free surface. The dimensionless parameters entering 
the above set of equations are the capillary number, the 
Marangoni number, the Galileo number, the Biot num- 
ber, and the Prandtl number: 
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and P = v jx- Here a is the surface tension, ctt = da/dT, 
g is the acceleration of gravity, q is the heat transfer rate, 
k is the thermal conductivity, x is the thermal diffusiv- 
ity, v and r\ are the kinematic and dynamics viscosity of 
liquid, respectively. 

Below we study evolution of a large-scale convection 
using the set of Eqs. ([T]) and @. 

Amplitude equations. We rescale the coordinates and 
the time as follows: 



X = ex, Y = ey, r = e t, 



(3) 



where e < 1 is the ratio of Hq to a typical horizontal 
lcngthscale. The temperature field is represented as T = 
-z + BP 1 +0(X,Y,t) + 0(e 2 ). 
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We assume large values of Ca and small values of Bi, 

Ca = e~ 2 C, Bi = e 2 p. (4) 

Thus we deal with the intermediate asymptotics between 
the conventional longwave mode, Bi = O (e 4 ) , § and the 
case of finite Bi |3| . These cases correspond to f3 = and 
P — » oo, respectively. 

Substituting the rescaled fields into Eqs. (UJ and © 
and applying the conventional technique of the lubrica- 
tion approximation (see 0), we arrive at 



h 3 

y vn 



V 

V • (W0) 

+v 



Mh 2 



V(6-h) 



Vj, 



(5) 



-(Vh) 2 -j3(6-h)+i-V(6-h) 



Mh 3 



V{9-h) 



(0) 



Here II = Gh — CV 2 h and V is a two-dimensional gradi- 
ent with respect to X and Y. 

Equations (|5|) and ([6]) form a closed set of the ampli- 
tude equations governing the nonlinear interaction of two 
well-known longwave modes: the Pearson's mode (h = 1) 
3 and the surface deformation-induced mode. (Note 
that the latter mode with 9 = const emerges only in the 
case of the conductive substrate Q.) Conductive state 
obviously corresponds to h = 9 = 1. 

Linear stability analysis. Substituting the perturbed 
fields /i = l+ £and6> = l + 9 into Eqs. © and ©, lin- 
earizing the equations for perturbations about the equi- 
librium, and representing the perturbation fields propor- 
tional to exp (At + ikX), one arrives at 
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FIG. 1: (a): Marginal stability curves M»(&) for G = 10: 
solid lines correspond to the monotonic mode, dashed ones 
- to the oscillatory mode; = 1, 10, 40 for lines 1, 2, and 
3, respectively, (b): The domain of oscillatory instability. 
The dashed vertical line marks the boundary of the longwave 
instability, Eq. ©. 



the limit C — is well studied [Tfj, for all computations 
we set C = 1 without loss of generality [ll| . The critical 
wavenumber materializing the minimum of the marginal 
stability curve, Eq. @, is 



(fcf»)) 



} \2 _ /3CG + ^/72(3CG (G + 72 - (iC) 
~ C (72 - (3C) ' 



For the oscillatory mode the marginal stability curve 
is determined by the expression 



M = 3 + G+ 



(10) 



The imaginary part of the growth rate for neutral per- 
turbations is 



A 2 + A 



a ,2 , G-M 

P + k 2 1 + 



+ y(/? + fc 2 )G- 




= 0, 



(7) 



where G = G + Ck 2 . Equation (JT)) possesses both real 
(monotonic instability) and complex (oscillatory instabil- 
ity) solutions. 

For the monotonic mode A = at the stability border, 
thus the marginal stability curve is given by 



M m = 



48 (P + k 2 ) G 



k 2 [72 + G 



(8) 



These marginal curves have a minimum at the finite val- 
ues of k only if 



PC < 72, 



(9) 



otherwise the minimal value, Mi m \ is achieved in the 
limit fc — > oo, i.e. the longwave mode is not critical. 
Hereafter we assume that the inequality (0) holds; since 



A, = Im(A) = —y (72 + G) (M m - M ), (11) 

i.e. the oscillatory mode is present only at M a (k) < 
M m {k). 

Minimization of the Marangoni number with respect 
to k gives 



Mj°) = 3 + G + 2^ZPC, fcW = ( M 



1/4 



(12) 



Examples of the marginal stability curves for these 
modes are shown in Fig.[TJa). Domains of monotonic and 
oscillatory instability are demonstrated in Fig.QJb). It is 
clear that the oscillatory mode is critical for PC > 17 A 
and G < 17.2. Take, for instance, a layer of water of 
thickness H = lO^cm. Then G w 0.1, Ca « 10 4 and 
Bi has to be approximately 10~ 3 in order to provide the 
required value of PC; this value seems achievable in ex- 
periments. 

Equations (fTU|) - (fT2")) indicate why the oscillatory mode 
has not been found earlier. As we have emphasized 



above, all previous studies deal with either G ^> 1 
or P = [|, or C = 
mode does not exist. 



In these cases the oscillatory 
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FIG. 2: (Color online). Pattern selection for the monotonic 
mode, (a) and (b) - the domains of stability for Rolls (marked 
with "R") and Squares ("S") on the square lattice. Solid 
(dashed) lines separate between supercritical and subcritical 
branching for Rolls (Squares). The latter domains are marked 
by "sub. R" ("sub. S"). Dotted lines separate domains of 
stability for Rolls and Squares. Dashed-dotted line in panel 
(a) is the locus of points N = 0; in the vicinity of this curve 
Eq. ()17[) holds. Diamond (circle) shows the threshold value 
Gi (G2) for pattern selection on the hexagonal lattice. 

Weakly nonlinear analysis. Monotonic mode. Here we 
study the nonlinear dynamics of perturbations at small 
supercriticality, M — Mc ~ 0, see Ref. @- To this end, 
we represent the primary part of the small perturbation 
of h in the form: 

71 

Z = ^2Aj exp (ikj ■ R) + c.c. (13) 

3=1 

where c.c. denotes complex conjugate terms and kj = 
k c m \ (The primary part of is expressed in terms of 
£.) The amplitudes Aj are functions of a slow time. 
For square (n — 2) and hexagonal (n = 3) lattices, the 
wave vectors are 

kx = fc c (l,0), k 2 =fc c (0,l) (14) 

andk! = fc c (l,0), k 2i3 = -fc c (-l,±V3), (15) 

respectively. 

For square lattice, the amplitude equations read 

Aj = (7 - K Q \A 3 \ 2 - KxSa) Aj, j = 1, 2, (16) 

where Sa — J2i l^'l 2 - Here the dot denotes the deriva- 
tives with respect to the slow time, and 7 ~ M — Mc 
is the real growth rate. The Landau constants, K n and 
K\ are real; they are cumbersome and thus are not pre- 
sented here. Results of the numerical calculations are 
shown in Fig. [2j One can readily see that supercriti- 
cal branching occurs only in two domains of parameters. 
These domains are situated either at rather small values 
of (3C, Fig. 0a), or at sufficiently small G, Fig. 02b). 
In the former case Rolls are selected everywhere except 
for a very small region shown in the inset. In the latter 
case Squares are selected everywhere excluding the small 
region where Rolls are stable. 



For hexagonal lattice, the resonant quadratic interac- 
tion results in the following amplitude equation: 

Ai = 1 A 1 - NA* 2 A* Z - (KolA^ 2 + K\Sa) A x , (17) 

and a similar equations for A 2t 3- (Hereafter the asterisk 
denotes the complex-conjugate terms.) Generally speak- 
ing, the quadratic term prevails over cubic ones, which 
leads to subcritical excitation of the hexagonal patterns 
through a transcritical bifurcation Q. However, N = 
at the dashed-dotted line shown in Fig. 02a) and in the 
vicinity of this line Eq. (IT71) becomes appropriate. 

Among the variety of possible patterns [7j, three are 
important. They are Rolls with A\ ^ 0, A2 = A3 = 
and two types of Hexagons with A\ = A2 = A3 = A: H + 
for A > and H~ in the opposite case. In the former 
case the flow is upward in the center of the convective 
cell, whereas in the latter case it is downward. 

Pattern selection on a hexagonal lattice is shown in 
Fig. 02a). At G < Gi s» 8.20 there are no stable solu- 
tions; the subcritical bifurcation occurs for Rolls and one 
branch of Hexagons (either H~ below or H + above the 
dashed-dotted line). At G\ < G < G 2 = 10 Rolls are still 
subcritical and unstable; stable Hexagons emerge only 
within the finite interval of supercriticality. Finally, at 
G > G2, H~(H + ) is stable within the interval of super- 
criticality, whereas Rolls become stable when M — Mj m ^ 
increases. 

To finalize the discussion of steady patterns, we briefly 
discuss the competition of patterns on the square and 
hexagonal lattices. It is clear that at the finite values of 
N , Hexagons emerge subcritically and no stable patterns 
can be found near the stability threshold. Therefore, 
weakly nonlinear analysis provides stable patterns only 
near the dashed-dotted curve shown in Fig. 02 a ) 1 where 
the competition between Hexagons and Rolls occurs. 

Weakly nonlinear analysis. Oscillatory mode. For the 
oscillatory mode the solution is presented in the form 

n 

S,=J2 (• l . ( ' k " + /<',< ' k ") e lAlT + c.c. (18) 
3=1 

Note that the pair (Aj , Bj ) corresponds to counter- 
propagating waves, which must be taken into account 
separately. The wavevectors for the square and hexago- 
nal lattices are given by Eqs. (IT41 and (|T3|). respectively. 

For square lattice, the equation governing the dynamics 
of the amplitudes Aj reads: 

Aj = [j-KolArf-K.lBtf -K 2 (S A + S B )]A 3 

-K 4 B*S AB ; j = 1,2, (19) 

where S B = Yli \ B i\ 2 > s ab = YJi A i B i- A similar pair 
of equations for Bj is obtained from Eqs. (fH)]) by replace- 
ment Aj -o- Bj. The Landau coefficients K\ (I = 0, 1, 2, 4) 
as well as the growth rate 7 are now complex- valued. 

Equations P^)l were studied in details in Ref. Q. Us- 
ing the results of that paper, we found that Traveling 
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FIG. 3: (Color online). Pattern selection for the oscillatory 
convection, (a) - square lattice: Domains of stability for TR 
(below the dashed line) and AR (to the left of the dotted 
line). Above the dashed line TR bifurcate subcritically. (b) 
- hexagonal lattice: Domains of stability for TR (below the 
dashed line and to the right of the dotted line) and TRa2 (be- 
tween the dotted and the solid line) are marked by "TR" and 
"TRa2", respectively. Above the dashed line TR bifurcate 
subcritically, to the left of the solid line TRa2 are subcritical. 



Rolls (TR), Ai 0, A2 = Bi,2 = can branch either su- 
percritically or subcritically [see Fig. EJa)], whereas the 
remaining patterns emerge through the direct Hopf bi- 
furcation; TR are selected in the domain of supercritical 
excitation. Alternating Rolls are stable within the small 
area marked by "AR" ; here depending on the initial con- 
dition the system either approaches AR or demonstrates 
the infinite growth of TR. 

For hexagonal lattice, the amplitude equation govern- 
ing the dynamics of the complex amplitudes Aj reads: 



Aj = [7 — if |A,-| 2 — Ki\Bj\ 2 — K 2 Sa — K 3 Sb] Aj 



-K4B* Sab, j — 1, 2, 3. 



(20) 



Three similar equations are obtained from Eqs. (|20[) by 
a replacement Aj f-> Bj. 

Analysis of the Hopf bifurcation for the above set of 
equations was performed in Ref. [s|, where eleven limit 
cycles were found and studied. Based on that paper, the 
results on pattern selection are presented in Fig. [HJb). 



The dashed line again separates direct and inverse Hopf 
bifurcations for TR, it is obviously the same as in the 
panel (a). However, for the hexagonal lattice, there ap- 
pears a competition between TR and Traveling Rectan- 
gles 2 (TRa2, A\ = B3 7^ 0, whereas all other ampli- 
tudes vanish). The latter pattern is stable in the domain 
marked by "TRa2" . The entire domain of supercritical 
bifurcation becomes smaller because TRa2 can bifurcate 
either supercritically or subcritically. 

Studying the competition between patterns on hexag- 
onal and square lattices, we found that the stability 
boundaries for both TR and TRa2 are the same as shown 
in Fig. GUb), whereas stability domain for AR nearly dis- 
appears. 

Conclusions. We studied the longwave Marangoni con- 
vection in a liquid layer heated from below; the heat flux 
at the substrate is specified. In such setup, an interaction 
of two well-known monotonic modes of longwave insta- 
bility, the Pearson's mode and the surface deformation- 
induced mode, can result in the emergence of a longwave 
oscillatory mode. However, the oscillatory mode has not 
been detected in spite of extensive numerical, analyti- 
cal, and experimental studies [l[ since the publication of 
Pearson's paper. We succeed in such analysis and point 
out the domain of parameters where the oscillatory mode 
exists, which can be reached in experiments. 

Moreover, we point out the domains of parameters 
where the convection emerges supercritically and hence 
either stationary or oscillatory terminal state with dis- 
torted surface is stable. This result is also very unusual, 
since only subcritical branching was found in the previ- 
ous studies [3, HI ■ 
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